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Abstract 

The paper presents the QCD description of the hard and semihard processes in the 
framework of the Wilson operator product expansion. The smooth transition between the 
cases of the soft and hard Pomerons is obtained. 
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The recent measurements of the deep-inelastic (DIS) structure function (SF) F2 by 
the HI [m and ZEUS collaborations open a new kinematical range to study proton 
structure. The new HERA data show the strong increase of F2 with decresing x. However, 
the data of the NMC ^ and i?665 collaboration ^ at small x and smaller is in the 
good agreement with the standard Pomeron or with the Donnachie-Landshoff picture 
where the Pomeron intercept: Op = 1.08, is very close to standard one. The interpritation 
of the fast changing of the intercept in the region of between = IGeV"^ and 

= lOGeV^ (see Fig. 3 in 0) is yet absent. There are the arguments in favour of that 
is one intercept (see 0) or the superposition of two different Pomeron trajectories, one 
having an intercept of 1.08 and the one of 1.5 (see Fig. 4 in |^). 

The aim of this article is the possible "solution" of this problem in the framework of 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation 0. It is good known (see, 
for example, M), that in the double-logarithmical approximation the DGLAP equation 



solution is the Bessel function, or exp y 0((5^)/n(l/x), where 0(Q^) is known Q^-dependent 
function^. However, we will seek the "solution"^ of DGLAP equation in the Regge form 
(we use the parton distributions (PD) multiplied by x and neglect the nonsinglet quark 
distribution at small x): 

fa{x,Q')^x-'Ux,Q^), {a = q,g) (a^ = 1 + <5) (1) 

where fa{x, Q"^) is nonsingular at x — > and /a(a;, Q^) ~ (1 — xY at s ^ iQ. The similar 
investigations were already done and the results are good known (see 0, [|1^-[|16])Q. The 
aim of this letter is to expand these results to the range where 5 ~ (and is not large) 



following to the observed early (see |]13|, |T^)0 method to replace the Mellin convolution 
by a simple product. Of course, we understand that the Regge behaviour (P is not in the 
agreement with the double-logarithmic solution, however the range, where 5 ~ and the 

values are nonlarge, is really the Regge regime and a "solution" of DGLAP equation in 
the form of (|1|) would be worthwhile. This "solution" may be understand as the solution 
of DGLAP equation together with the condition of its Regge asymptotic at a; — 0. 

Consider DGLAP equation and apply the method from to the Mellin convolution 
in its r.h.s. (in contrast with standard case, we use below a(Q^) = as((5^)/(47r)): 

^/a(3:,t) = - ^ X! 7ai(a,a;) (g) /a(x,t) {a,b) = {q,g) 

i=a,b 

= - ^ E 7a^(«, 1 + 5)fa{x,t) + 0{x'-') {^ab{a,n) = a7SV) + cy'lH^ (n) + , 

i=a,b 



^More correctly, (p is Q^-dependent for the solution of DGLAP equation with the boundary condition: 
fa{x, Qq) — Const at X ^ 0. In the case of the boundary condition: fa{x, Qq) ~ exp ^Jln(\/x), (j) is lost 
(see [||) its Q^-dependence 

■^We use the terrain "solution" because we will work in the leading twist approximation in the range of 
Q^: > IGeV^, where the higher twist terms may give the sizeable contribution (see, for example, pO|). 
Moreover, our "solution" is the Regge asymptotic with unknown parameters rather then the solution of 
DGLAP equation. The parameters are found from the agreement of the r.h.s. and l.h.s. of the equation. 

^Consideration of the more complicate behaviour in the form x^^ {ln{l/x))^l2g{^y <j)ln{l/x)) is given 
in 0] and will be considered in this content in the forthcomming article pj} ] 

^n the double-logarithmical approximation the similar results were obtained in [T^ 

^The method is based on the earlier results 
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where t = ln{Q'^/h?'). The 7ab(a,x) are the sphting functions corresponding to the 
anomalous dimensions (AD) 7a6(a, n) = (ixx"'~^7ab(a, x). Here the functions 7afe(a, 1 + 
6) are the AD 7ab(a,?2) expanded from the integer argument "n" to the noninteger one 
"1 + 6" . The functions 7a6(«, 1 + 5) (marked lower as AD, too) can be obtained from the 
functions 7afe(a, 1 + 5) replacing the term 1/S by the one 1/6: 



1 

6 



1 

6 



t(i - v{x,S)x' 



(3) 



This replacement (^ is appeared very naturally from the consideration the Mellin con- 

Q) and preserves the smooth and nonsingular transition to the 



volution at X ^ 
case 5 = 0, where 



see 



1 



In g(x) 

X 



(4) 



The concrete form of the functions (p{x, 6) and g{x) depends strongly on the type of 
the behaviour of the PD fa{x, Q^) at x and in the case of the Regge regime (|lD they 
are (see P, |I1): 



¥^(^,5) = ^r/i^y^\/^ and gix) = + (5) 

where r(z/ + 1) and \l/(i^ + 1) are the Eulerian F- and \l/-functions, respectively. As it can 
be seen, there is the correlation with the PD behaviour at large x. 

If 6 is not small (i.e. x^ » 1), we can replace 1/6 to 1/6 in the r.h.s. of Eq.(^ and 
obtain its solution in the form (hereafter to = ti^Q"^ = Ql)): 



fa{x,t) ^ Ma{l + 6,t) 

faix,to) M,(l + 5,to)' ^' 

where M„(l+5, t) is the analytical expansion of the PD moments Ma{n, t) = dxx^^^^ fa{x, t) 
to the noninteger value "tt, = 1 + 5" . 

This solution is good known one (see [0| for the first two orders of the perturbation 



theory, ||T5[ for the first three orders and containing a resummation of all orders 



respectively) . Note that recently the fit of HERA data was done in with the formula 
for PD fq{x,t) very close [] to (|^) and the very well agreement (the P^r degree of 
freedom is 0.85) is found at 5 = 0.40 ± 0.03. There are also the fits of the another 
group using equations which are similar to (|^) in the LO approximation. 

The news in our investigations are in the follows. Note that the Q^-evolution of Ma{l + 
6,t) contains the two: "+" and "-" components, i.e. M„(l + 6,t) = Ei=±M*(l + 5,t), 



^The used formula (Eq.(2) from |T1]) coincides with (^) in the leading order (LO) approximation, if 
we save only fg{x, Q^) in the r.h.s. of (|) (or put jqq = and jqg = formally). Eq.(||) and Eq.(||) from 
[jiol have some differences in the next-to- leading order (NLO), which are not very important because they 
are corrections to the a-correction. 
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and in principle the every component evolves separately and may have the independent 
(and not equal) intercept. Here for the simplicity we restricte ourselves to the LO analysis 
and give NLO formulae lower without large intermediate equations. 

1. Consider DGLAP equation for the "+" and "— " parts (hereafter s = ln{lnt/lnto)): 



0(x 



l-5\ 



(7) 



where 



7± 



(733 + 7gg) ± \l{lgg-lqq) + ^Iqglgq 



are the AD of the "±" components (see, for example, ||21|| ) 

The "— " component 7- (a, 1 + 5 J) does not contain the singular term (see [|13|, and 
lower) and its solution have the form: 



- e-^-(^-^^-^ where d + 



fai^M ' ^ 2(3, 

The "+" component 7+(a, 1 + 5+) contains the singular term and f^{x,t) have the 
solution similar (|^) only for >> 1: 

= e"'^+(l+^+)^ if » 1 (9) 

The both intersepts 1 + 5+ and 1 + 5_ are unknown and should be found, in principle, 
from the analysis of the experimental data. However there is the another way. From the 
small (and small x) data of the NMC and EQQb collaboration 0] we can conclude 
that the SF F2 and hence the PD faix-iQ"^) have the flat asymptotics for x ^ and 

~ (1 2)GeV'^. Thus we know that the values of 5+ and 5_ is approximately zero at 

Consider Eqs.(0) with 5± = and with the boundary condition fa{.x,Ql) = Aa at 
Qq = IGeV"^. For the "— " component we already have the solution: the Eq.(|^) with 
5_ = and (i_(l) = 16//(27/3o), where / is the number of the active quarks and are 
the coefficients in the a-expansion of QCD /3-function. For its "+" component Eq.(|^) can 
be rewritten in the form (hereafter the index 1 + 5 will be omitted in the case 5^0): 



(10) 



where 7+ and 7_,_ are the coefficients of the singular and regular parts at 5 ^ of AD 
7+(l + 5): 

7+(l + 5) = ^+1 + ^+ = -24, 7+ = 22 + ^ 

The solution of Eq.(lO) is 
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/+(a;,t) = Atx^^^e-~'^\ (11) 

where 

7+ 4 - 1 /_ \ 4 101 

" 2^ " -3' ^+ " WS^^ " ^"'^'^^ " 3^^"^ + ^ 
Herefter the symbol ~ marks the case / = 3. 

As it can be seen from (^) the flat form 5+ = of the "+" -component of PD is very 
nonstable from the (perturbative) viewpoint, because d{5+)/ds ^ 0, and for Q"^ > Ql 
we have already the nonzero power of x (i.e. pomeron intercept Up > 1). This is in the 
agreement with the experimental data. Let us note that the power of x is positive for 
Q2 < Ql that is in principle also supported by the NMC data, but the use of this 
analysis to < IGeV^ is open the question. 

Thus, we have the DGLAP equation solution for the "+" component at is close to 
Qq = IGeV"^, where Pomeron starts in its movement to the subcritical (or Lipatov E^, |23|) 



regime and also for the large Q^, where pomeron have the Q^-independent intercept. In 
principle, the general solution of (|^) should contain the smooth transition between these 
pictures but this solution is absent |. We introduce the some "critical" value of Q^: Ql, 
where the solution (|^) is replaced by the solution (^). The exact value of Ql may be 
obtained from the flt of experimental data. Thus, we have in the LO of the perturbation 
theory: 



where 



Ux,t) = f-{x,t) + f+{x,t) 
fai^^t) = A~ exp{-d^s) 

^ , Atx^-^-^expi-d+s), ifQ'<Ql 

' ^ ^ /+(x,gexp(-rf4i + ^c)(.-.c)), ifg^>g^ ^ ' 



tc 

and A 



t{Qi), = s{Qi) 
(1 - a)Ag + aAg, 
= A. - At 



A'„ 



aAg — eAq 



(13) 



and the values of the coefficients a, a and e may be found, for example, in |^ 
Using the concrete AD values at 5 = and / = 3, we have 



A-1 



1 



AAa + 9 A, 



27/n(i)-,(z.)-^ 



^9 + - 



108 

4 9A„ - A 



27/n(i) 



(14) 



^The form exp(^-s7+(l + (5)/(2/3o)) coincides with the both solution: Eq.(g) if x'^+ » 1 and Eq.(|ri] 
when ^ = but it is not the solution of DGLAP equation. 



4 



Thus, the value of the "+" component of the quark PD is suppressed logarithmically 
that is in the qualitative agreement with the HERA parametrizations of SF F2 (see 
H, 13) (in the LO F2{x,Q^) = (2/9) fg{x,Q^) for / = 3), where the magnitude 



connected with the factor x'^ is 5 -^ 10% from the flat (for x — > 0) magnitude. 

2. By analogy with the subsection 1 and knowing the NLO Q^-dependence of PD 
moments, we obtain the following equations for the NLO Q^-evolution of the both: "+" 
and "— " PD components (hereafter s = ln{a{Ql)/a{Q'^)),p = a{Ql) — a(Q^)): 

faix,t)= A- exp (-rf_S - 

r+(^ ,^ ^ / Atx^^-'^^^UP) exp (-rf+S - dl^p), if g2 < Q2 

^"^'^ \ /+(x,gexp(-rf+(i + 5,)(s-5,)-<+(i + 5,)(p-p,)), ifg2>gr ' 

where 



siQl), Pc = p{Ql), tto 



«(Qo) 



1 - aoKl)A^ + a^K^A"^ 



d\Aln{-) - q{u)) + d 



X 



and Kl 



7±: 



Ki 



Kl 



_ 7±± 

(0) 

7t - 7w 



(16) 



2/3o + 7±-7t' 

The NLO AD of the "±" components are connected with the NLO AD 7^),^ The 



corresponding formulae can be found in |21 



Using the concrete values of the LO and NLO AD at 5 = and / = 3, we obtain the 
following values for the NLO components from (p^) , (p^ (note that we remail only the 
terms ~ 0(1) in the NLO terms) 



and 



dl 



d^ 



16 
81 



2C(3) + 9C(2) 



779 



108 



1.97, dL 



d'i__ + ^ 2.32 



2800 -q 

81 ' +■ 
1180 , 



81 



32 



d. 



C(3) + |!C(2) 



372607 



69984 



-67.82 



953 



- 12C(2) ?^ -52.26 



(17) 



20 

y 



^9 + g A 



+ 



1 4Ag(l - 7.67ao) + 9Ag{\ - 8.71ao) 



27 



/n(i) - q{v) - 



85 
108 



and A„ 



Aa — A 



9Ag - A^ 



27/n(i)-f?(z/)- 



85 
108 



/ 692 , 
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It is useful to change in Eqs.(|T5|)-(|TB|) from the quark PD to the SF F2(x, Q^), which 
is connected in NLO approximation with the PD by the following way (see pT|): 



l + a{Q')B,{l + 6)]6if,{x,Q') + a{Q')B,{l + 6)6'Jg{x,Q' 



(19) 



where 5^ = J2i=i If =< ^} > is the average charge square of the active quarks: 5^ = (2/9 
and 5/18) for / = (3 and 4), respectively. The NLO corrections lead to the appearence in 
the r.h.s. of Eqs.(|l5D of the additional terms ^1 + aB-j-^ /(''"''' '^oB±j and the necessarity 

to transform to = F^{x, Q^) into the input parts. The final results for F2{x, Q^) 
are in the form: 



F2(x,t) 

F+ix,t) 
where 



F^-{x,t) + F+{x,t) 

C- exp {-dSs - dt_p){l + aB-)l{l + «o5") 
f exp {-d+~s - rf++p)(l + aB+)/{l + ao5+) 

F}{x,t^) exp (-rf+(l + 5c) (5 - Sc) - c?++(l + 4)(p - P, 

\ + aB+{l + 5,)]/(l + a,B+{l + 4) 



B^ = B, + ^Bg, C± = A^H + a^B' 



with the substitution of Aq hy C = F2{x, Ql) into Eq.(|18]) according 



if Q' < Ql 

(20) 

if > Ql 



C = (l + aoBq)6lAq + aoBgd^Ag, 



(21) 



For the gluon PD the situation is more simple: in Eq.(|IB]) it is necessary to replace 
Ag by C according (pi]). 

For the concrete values of the LO and NLO AD at 5 = and / = 3, we have for 
Q^-evolution of F2{x, Q"^) and the gluon PD: 



F2{x,t) 
F^-{x,t) 

F^{x,t) 



= F^-{x,t) + F^{x,t), fg{x,t) = f-{x,t) + f^{x,t) 
= C- exp(-^5-1.97p)(l-^a)/(l-^ao) 



= < 



^+^(_|s+2|20p) (_|(^,(z.) + mys + C^q{v) - 67.82)p 
(l + 6[/n(i) - q{u) - ii]a)/(l + 6[Zn(i) - q{u) - i|]«o' 
F}{x, Q exp (-c/+(l + 6,){s - Sc) - rf^+(l + S,){p - p^^ 

'l + a5+(l + 4)) / (l + «c5+(l + 4)) , 
-|5-2.32p)(l-^a)/(l-^ao) 



if < 

(22) 



if > g: 



exp 
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where 



' exp {-l{g{u) + + C-^g{u) - 52.26)p 

(l + 6[/n(i) - giu) - l]a)/(l + 6[/n(i) - giu) - l§]ao), if Q' < Q 
f+{x, Q exp + - Se) + + - Pc)] 

1 + aB+{l + 4)) / f 1 + + 5,)) , if Q2 > q2 



23) 




v4^ + 2C 



+ 



Agil - 9.74«o) + 2^(1 - 7.82ao) 



85 
108 



2C 



2 2A,(l + f^«o)-C^(l + ^ao 



27 



and A„ 



A, 



(24) 



(25) 



Let us give some conclusions following from Eqs.(l2^)- (|25|) . It is clearly seen that 
the NLO corrections reduce the LO contributions. Indeed, the value of the subcritical 
Pomeron intercept, which increases as ln{ao/(y) in the LO, obtaines the additional term 
~ (ao ~ cfc) with the large (and opposite in sign to the LO term) numerical coefficient. 
Note that this coefficient is different for the quark and gluon PD, that is in the agreement 
with the recent MRS{G) fit in |26| and the data analysis by ZEUS group (see [Q). The 
intercept of the gluon PD is larger then the quark PD one (see also [EB, P3|). However, the 



effective reduction of the quark PD is smaller (that is in the agreement with W.-K. Tung 
analysis in [^), because the quark PD part increasing at small x obtains the additional 
(~ ao but not ~ 1/lnx) term, which is important at very small x. 

Note that there is the fourth quark threshold at Q'^f^^ ~ lOGeV^ and the Q^/j value may 
be larger or smaller to Ql one. Then, either the solution in the r.h.s. of Eqs. (p0|j2^j23|) 



3 to ones at / = 4. The 
(see also the recent experimental test of the flavour 



before the critical point Qc and the one for Q > contain the threshold transition 
where the values of all variables are changed from ones at / 
aiQ'^) is smooth because — > A^ti 
independence of strong interactions into 

For simplicity here we suppose that Qff^ = Ql and all changes initiated by threshold 
are done authomatically: the flrst (at Q"^ < Ql) solutions contain / = 3 and second (at 
Q'^ > Ql) ones have f = 4:, respectively. For the "— " component we should use Qff^ = Q^, 
too. 

Note only that the Pomeron intercept = 1 — {d+s + d\^p) increases at Q^ 
because 



Q 



2 

thy 



lAASiQl, Ql) 



38AlpiQl,Ql), 



81 



if Q' < Ql 
if Q^ > Ql 



that agrees with results [29 obtained in the framework of dual parton model. The differ- 



ence 



Aap = 0.115(QLQo)-3.55p(QLQo) 
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dependes from the values of and Q^. For Qff^ = lOGeV^ and Ql = IGeV"^ it is very 
small: 

Aap = 0.012 

3. Let us resume the obtained results. We have got the DGLAP equation "solution" 
having the Regge form (|1]) for the two cases: at small {Q"^ ~ IGeV"^), where SF and 
PD have the flat behaviour at small x, and at large Q^, where SF F2{x, Q^) fastly increases 
when X ^ 0. The behaviour in the fiat case is nonstable with the perturbative viewpoint 
because it leads to the production of the subcritical value of pomeron intercept at larger 
and the its increase (like 4/3 ln{a{Ql) / a^Q"^) in LO) when the value increases]^. 
The solution in the Lipatov Pomeron case corresponds to the well-known results (see 
T3| , ITsj |19[] ) with Q^-independent Pomeron intercept. The general "solution" should 



contains the smooth transition between these pictures. Unfortunately, it is impossible 
to obtain it in the case of the simple approximation (|I]), because the r.h.s. of DGLAP 
equation (|^) contains the both: ~ x~^ and ~ Const, terms. As a result, we used two 
above "solutions" gluing in some point Q^. 

Note that our "solution" is some generation (or a application) of the solution of 
DGLAP equation in the momentum space. The last one have two: "+" and "— " com- 
ponents. The above our conclusions are related to the "+" component, which is the 
basic Regge asymptotic. The Pomeron intercept corresponding to "— " component, is Q^- 
independent and this component is the subasymptotical one at large Q^. However, the 
magnitude of the "+" is suppressed like l/ln{l/x) and a{Ql), and the subasymptotical 



"— " component may be important. Indeed, it is observed experimentally (see |2J, psp. 
Note, however, that the suppression ~ o:{Qo) is really very slight if we choose a small 
value of Qo- 



Our "solution" in the form of Eqs.([2^)-(P3|) is in the very well agreement with the 



recent MRS{G) fit [||] and with the results of |T9l at = i^GeV^. As it can be seen 
from Eqs.(p^),(p3l), in our formulae there is the dependence on the PD behaviour at large 
X. Following to ||32| we choose u = 5 that agrees in the gluon case with the quark counting 
rule This u value is also close to the values obtained by CCFR group |Q {v = 4) 



and in the last MRS{G) analysis {u = 6). Note that this dependence is strongly 
reduced for the gluon PD in the form 

f,ix,Ql) = A,(z/)(l-xr, 

if we suppose that the proton's momentum is carred by gluon, is zz-independent. We used 
Ag{5) = 2.1 and F2(x, Ql) = 0.3 when x ^ 0. 

For the quark PD the choise z/ = 3 is more preferable, however the use of two different 
u values complicates the analysis. Because the quark contribution to the "-I-" component 
is not large, we put z/ = 5 to both: quark and gluon cases. Note also that the variable 
z/((5^) have (see |^|) the Q^-dependence determinated by the LO AD 7^5. However this 



Q^-dependence is proportional s and it is not important in our analysis. 

Starting from Ql = IGeV"^ (by analogy with |^) and from Ql = 2GeV^, and using 



two values of QCD parameter A: more standard one (A:^ = 200 MeV) and (A^^= 



^The Pomeron intercept value increasing with was obtained also in pC 
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255 MeV) obtained in 
"intercepts" 6a = — {d 
if g2 = 1 GeV^ 



we have the following values of the quark and gluon PD 
s + d\^a) (here is marked as A): 



if Ql 















A = 200MeF 


A = 200MeV 


A = 255MeV 


A = 255MeV 


4 


0.191 


0.389 


0.165 


0.447 


10 


0.318 


0.583 


0.295 


0.659 


15 


0.367 


0.652 


0.345 


0.734 


2 GeV^ 


Q' 


W) 




W) 






A = 200MeV 


A = 200MeV 


A = 255MeV 


A = 255MeV 


4 


0.099 


0.175 


0.097 


0.198 


10 


0.226 


0.368 


0.227 


0.410 


15 


0.275 


0.438 


0.278 


0.486 



Note that these values of 6a are above the ones from p6|. Because we have the second 
(subasymptotical) part, the effective our "intercepts" have the smaller values. 

As a conclusion, we note that BFKL equation (and thus the value of Lipatov Pomeron 
intercept) was obtained in in the framework of perturbative QCD. The large-Q^ 
HERA experimental data are in the good agreement with Lipatov's trajectory and thus 
with perturbative QCD. The small data agrees with the standard Pomeron intercept 
ftp = 1 or with Donnachie-Landshoff pisture: ctp = 1.08. Perhaps, this range requires 
already the knowledge of nonperturbative QCD dynamics and perturbative solutions 
(including BFKL one) should be not applied here directly and are corrected by some 
nonperturbative contributions. 

In our analysis Eq.(|I]) can be considered as the nonperturbative (Regge-type) input 
at Qq ~ IGeV"^. Above Qq the PD behaviour obeys DGLAP equation, Pomeron moves 
to the subcritical regime and tends to its perturbative value. After some Ql, where its 
perturbative value was already attained, Pomeron intercept saves the permanent value. 
The application of this approach to analyse small x data invites futher investigation. 
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